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I reanalyze the hydrodynamic theory of fluid, polar ordered flocks. I find new linear terms in the 
hydrodynamic equations which slightly modify the anisotropy, but not the scaling, of the damping of 
sound modes. I also find that the nonlinearities allowed in equilibrium do not stabilize long ranged 
order in spatial dimensions d = 2; in accord with the Mermin- Wagner theorem. Nonequilibrium 
nonlinearities do stabilize long ranged order in d = 2, as argued by earlier work. Some of these were 
missed by earlier work; it is unclear whether or not they change the scaling exponents in d = 2. 

PACS numbers: 05.65.+b, 64.70.cij, 87.18.Gh 



Flocking [l[ - the coherent motion of large numbers 
of self-propelled entities, - spans a wide range of length 
scales: from kilometers (herds of wildebeest) to microns 
(microorganisms 0, H ; mobile macromolecules in living 
cells 0,11). It is also [6| a dynamical version of ferroma_; 
netic ordering. A "hydrodynamic" theory of flocking 
shows that, unlike equilibrium ferromagnets Q, flocks 
can spontaneously break a continuous symmetry (rota- 
tion invariance) by developing long-ranged order, (i.e., a 
non-zero average velocity {v{r,t)) ^ 0) in spatial dimen- 
sions d = 2, even with only short ranged interactions. 

The mechanism for this apparent violation of the 
"Mermin- Wagner" theorem Q is fundamentally nonlin- 
ear. A number of nonlinear terms in the hydrodynamic 
equations of motion become "relevant" , in the renormal- 
ization group (RG) sense, as the spatial dimension d is 
lowered below 4, leading to a breakdown of linearized hy- 
drodynamics which suppresses fluctuations enough to 
stabilize long-ranged order possible in d = 2. 

In this paper, I revisit the formulation of the hydrody- 
namic theory of what I'll call "fluid, polar ordered" flocks, 
by which I mean flocks that are spatially homogeneous, 
on average, and havc_(u(r, t)) ^ 0. I find a few differences 
with the results of [7|. Some of these are minor: a few 
linear terms, that produce only minor modifications of 
the damping of the propagating sound modes predicted 
in Q, were missed in that earlier work. 

My more important conclusions concern the scaling 
laws of two dimensional flocks. It was originally argued 
0] that the exponents characterizing the scaling of fluc- 
tuations in flocks that results from the breakdown of hy- 
drodynamics could be determined exactly in d = 2. In 
this paper, I will argue that those arguments were incor- 
rect, because they neglected certain other, equally im- 
portant, symmetry- allowed nonlinearities in the hydro- 
dynamic equations. These additional nonlinearities in- 
validate the earlier arguments, and render it impossible 
to determine the exact scaling laws in d = 2, or, indeed, 
in any spatial dimension d < 4. 

// these new nonlinearities should prove to be irrele- 



vant, in the RG sense, in d = 2, then the exact expo- 
nents predicted by would, in fact, hold in d = 2. At 
the moment, however, there is no compelling theoretical 
argument that they are irrelevant, though there is also 
none that they are not. 

All of these non-linearities involve density fluctuations. 
Hence, in systems in which density fluctuations are sup- 
pressed, it is possible to obtain exact exponents in d = 2. 
One class of such systems - flocks with birth and death 
- has been treated elsewhere (loj: others, such as in- 
compressibl e sy stems and systems with long-ranged 
interactions [HI, will be addressed in future work. [l3| 

The new treatment presented here correctly predicts 
that one naively relevant nonlinearity in the flocking hy- 
drodynamic equations that is allowed even in equilibrium 
systems [l3l does not lead to any corrections to scal- 
ing (or, indeed, to any qualitatively new long-wavelength 
physics whatsoever); this means that the equilibrium sys- 
tems described by such a model does not exhibit long- 
ranged order in d = 2 (in accord with the Mermin- 
Wagner theorem Q). 

My discussion here is limited to "ordered" flocks mov- 
ing on a substrate: i.e., one in which the flocking or- 
ganisms spontaneously pick a direction to move together 
via purely short-ranged interactions that make neighbors 
tend to follow each other, but which do not pick out 
any a priori preferred direction for this motion. That is, 
the flocking spontaneously breaks rotation invariance, as 
equilibrium ferromagnetism does. Flocks moving with- 
out a substrate conserve momentum, and so have a 
very different hydrodynamics, which has been considered 
elsewhere [IB}: I will not discuss these here. One specific 
realization of a flock on a substrate is the Vicsek algo- 
rithm Qin its ordered state. 

The hydrodynamic theory describes the flock by con- 
tinuous, coarse grained number density p(r, t) and veloc- 
ity v{r,t) fields. The hydrodynamic equations of motion 
governing these fields can in the long-wavelength limit 
can be written down purely on symmetry grounds 0], 
and are: 
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dtv + Xi{v-V)v + X2{V ■v)v + X3^{\v\^) = av- l3\v\^v-^Pi -v(^v-VP2^ + DbV {V ■ v) + Dt^^v + D2{v ■ Vfv + f 



(1) 



dtP + V- (vp) = (2) 

where all of the parameters Xi{i = 1 3), a, /?, 
Db.t.2 and the "isotropic Pressure" P{p, \v\) and the 
"anisotropic Pressure" P2(p, l^l) are, in general, functions 
of the density p and the magnitude \v\ of the local veloc- 
ity. It is useful to Taylor expand Pi. 2 and P2 around the 
equilibrium density po: 

Here /3, Db-, -D2 and Dt are all positive, and a < in 
the disordered phase and a > in the ordered state (in 
mean field theory). 

The a and /? terms simply make the local v have a 
nonzero magnitude vq = in the ordered phase, 

where a > 0. Db,t,2 are the diffusion constants (or vis- 
cosities) reflecting the tendency of a localized fluctuation 
in the velocities to spread out because of the coupling 
between neighboring "birds" . The / term is a random 
driving force representing the noise. It is assumed to be 
Gaussian with white noise correlations: 

< Mr, t)fj{r\ t') A%<5'^(f - ?)5{t - t') (4) 

where A is a constant, and i , j denote Cartesian compo- 
nents. The pressure P tends, as in an equilibrium fluid, 
to maintain the local number density p{r) at its mean 
value Po, and 5p = p — po- The "anisotropic pressure" 
P2(p, \v\) in equation ([ij is only allowed due to the non- 
equilibrium nature of the flock; in an equilibrium fluid 



such a term is forbidden, since Pascal's Law ensures that 
pressure is isotropic. In the nonequilibrium steady state 
of a flock, no such constraint applies. In earlier work 0, 
this term was ignored. Here I will show that this term 
changes none of the predictions of the hydrodynamic the- 
ory. 

The flnal equation ([2]) is just conservation of bird num- 
ber: we don't allow our birds to reproduce or die on the 
wing. The interesting and novel results that arise when 
this constraint is relaxed by allowing birth and death 
while the flock is moving will be discussed elsewhere [lol|. 

The hydrodynamic model embodied in equations ([!]). 
([2]), dS]), and (g]) is equally valid in both the "disor- 
dered" (i.e., non- moving) (a < 0) and "ferromagncti- 
cally ordered" (i.e., moving) (a > 0) state . Here I 
am interested in the "ferromagnetically ordered" , broken- 
symmetry phase which occurs for a > 0. In this state, 
the velocity field can be written as: 



V = voX\\ + Sv — {vo + Sv\\)x\\ + v± , (5) 

where voX\\ =< v > is the spontaneous average value of 
V in the ordered phase, and the fluctuations Sv\\ and v± 
of V about this mean velocity along and perpendicular 
to the direction of the mean velocity arc assumed to be 
small. Indeed, I will be shortly be expanding the equation 
of motion ([T]) in these quantities. Taking vq — as 

discussed above [HI , and taking the dot product of both 
sides of equation ([T]) with v itself, I obtain: 



- (dt\vf + (Ai +2A3)(i' • V)|z;p) + A2(V ■ v)\vf = (a - l3\vf)\vf -v-VP- \v\^v ■ V P2 + Dbv ■ y{V ■ v) 

+ Dtv-V^v + D2V- (^{v-Vfv^ +v- f . (6) 
I 



In this hydrodynamic approach, we are interested only in 
fluctuations dv{r, t) and 6p{f, t) that vary slowly in space 
and time. (Indeed, the hydrodynamic equations ^ and 
^ are only valid in this limit). Hence, terms involv- 
ing space and time derivatives of 5v{f,t) and Sp{r,t) are 



always negligible, in the hydrodynamic limit, compared 
to terms involving the same number of powers of fields 
without any time or space derivatives. 

Furthermore, the fluctuations 6v{'r,t) and 5p{f,t) can 
themselves be shown to be small in the long-wavelength 
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limit. Hence, we need only keep terms in equation ([6]) up 
to linear order in Sv{r, t) and Sp{f, t). The v- f term can 
likewise be dropped, since it only leads to a term of order 
v±f\\ in the vj_ equation of motion, which is negligible 

(since v± is small) relative to the f± term already there. 

These observations can be used to eliminate many of 
the terms in equation (|6l), and solve for the quantity 

U^iaip,\v\)~l3ip,\v\)\v\^) ; (7) 

I 



This can be made into an equation of motion for v± in- 
volving only v±{r,t) and 5p(f,t) by projecting perpen- 
dicular to the direction of mean flock motion x\\, and 
eliminating 5v\\ using equation® and the expansion 

U « -Ti (^5v\\ + - T^Sp , (10) 

where I've defined 

I 



where I've kept only linear terms on the right hand side 
of this equation, since the non-linear terms are at least 
of order derivatives of \v±\'^ , and hence negligible, in the 
hydrodynamic limit, relative to the \v±\'^ term explicitly 
displayed on the left-hand side. 

This equation can be solved iteratively for 5v\\ in terms 

I 



the solution is: 

U = \2V ■V + V-VP2 + —d\\5p+^ {dt+i2d\0 ^f||,(8) 

Vq ZVq 

where I've defined 72 = (Ai + 2A3)uo- Inserting this ex- 
pression (lU for U back into equation ([6]) (where U ap- 
pears by virtue of its definition ([7]), I find that P2 and A2 
cancel out of the v equation of motion, leaving 



(9) 



with, here and hereafter , super- or sub-scripts denoting 
functions of p and \v\ evaluated at p = po and \v\ = vq. 
I've also used the expansion ([5]) for the velocity in terms 
of the fluctuations 5v\\ and v±_ to write 

\v\=vo + Sv\\+^-^ + 0{dvl\v^\') , (12) 

and kept only terms that an RG analysis shows to be 
relevant in the long wavelength limit. Inserting (|10p into 
(HI) gives: 



(13) 

I 

of v±, 5p, and its derivatives. To lowest (zeroth) order 
in derivatives, ~ —j^Sp. Inserting this approximate 
expression for into equation (|13p everywhere i5w|| ap- 
pears on the right hand side of that equation gives (jDy 
to flrst order in derivatives: 



dtv + Xi(v-V)v + AsVdwP) = ^v(d»6p) -VP + DbV(V ■ v) + DyV^w + L»2(w • V)^^ 

Vq 



2vo 



{dt + 729||) (5u|| 



v + f 



Ti Sv 



'II 



\v±\- 
2vo 



r2(5p = A2V_L • V±_ + A29||(5l'|| + 



(0-2, ifo + (^1,1) 



9\\ Sp + TT- {^t + l2d\\ ) 5v\\ 
II 2vr 11/11 



~ TTOtSp + — - — V_L • Vi_ - — 



where A4 is a constant related to the constants in equation Inserting ([5]), (jl2[) . and (jl4p into the equation of mo- 
([12]). tion ^ for v, and projecting that equation perpendicular 
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to the mean direction of flock motion i|| gives, neglecting "irrelevant" terms: 
I 



+ £>BeffV± fVi • + D^Vlv4_ + D?,df^V4_ + lytA^ ±Sp + unduV ±Sp + (15) 



r 



where Dg^g,D'^^, Si,2,3,4, Cq, and z/^ || are all constants Using ([5]) and in the equation of motion ^ for p 
expressible in terms of the parameters in equation (HH). gives, again neglecting irrelevant terms: 



dtSp + Po'^j 



{vj_5p) + V2d\\6p ■ 



Dp\\d'^^6p + Dpyd\\ 



-iTi) +w^dtd\\5p + w2d\\{5p') + ^d\\{\v^\') , 



2vo 



(16) 



where W2, wi,2, and -Dp||,p„ are also all constants express- 
ible in terms of the parameters in equation (|14p . 

I will henceforth focus my attention on the uniform 
ordered state, in which all of the diffusion constants Dp\\ , 
Dpv, D^BcS^ ^^'i positive. 

The linearized forms of equations ([T5|) and (|16P differ 
from the corresponding equations considered in [7| only 
in the Ut \\ terms in equation (|15p . and the -Dp || and Dp^ 
terms in equation (|16p . These prove [13[ to lead only to 
minor changes in the propagation direction dependence, 
but not the scaling with wavelength, of the damping of 
propagating sound modes predicted in Q ; their direction- 
dependent speeds are unaffected. 

The non-linear terms in equations (|15p and (|16p are 
more significant. Several of these, specifically all of those 
involving || and t derivatives, were missed by 0, in part 
because potential density dependences of various param- 
eters were missed, and in part because of subtle mistakes 
in eliminating the fluctuations 5v\^\ of the velocity along 
the mean direction of motion. 

Equally noteworthy are the non-linear terms that are 
missing from equations (jl5|) and (jl6p : all nonlinearitics 
arising from the anisotropic pressure P2 and the A2 non- 
linearity. This in particular has the very important con- 
sequence of saving the Mermin- Wagner theorem. This 
is because the A2 term is allowed even in equilibrium 
systems The incorrect treatment in suggested 

that this term by itself could stabilize long-range order 
in d = 2. Given that this term is allowed in equilib- 
rium, this would imply that the Mermin- Wagner theorem 
would fail for such an equilibrium system. The correct 
treatment I've done here shows that this is not the case: 
the A2 term by itself cannot stabilize long ranged order in 



d = 2, since the non-linearities associated with it drop out 
of the long-wavelength description of the ordered phase. 

Returning now to the non-linearities in equations (jisp 
and (flG]) that were missed by , I note that all of them 
become relevant, in the renormalization group (RG) 
sense (l7j. for spatial dimensions d < dc, (where dc is 
the critical dimension below which nonlinear effects be- 
come important; here dc = 4), according to simple power 
counting. This invalidates the arguments given in Q 
that gave exact exponents in d = 2, since those argu- 
ments were based, inter alia, on the assertion that all of 
the relevant nonlinearities could, in d = 2, be written 
as total 1- derivatives. This assertion is clearly not true 
of the terms involving || and t derivatives in equations 
([T5|) and (|T6| : hence, those terms completely invalidate 
the arguments leading to the exact exponents. I have 
been unable to come up with alternative arguments that 
give exact exponents in the presence of these additional 
terms. 

Now, if these additional nonlinearities were irrelevant 
in d = 2 under a full dynamical RG, then the exact ex- 
ponents of would be correct in d = 2. 

There is a precedent for this (that is, for terms that 
appear relevant by simple power counting below some 
critical dimension dc actually proving to be irrelevant 
once "graphical corrections" -i.e., nonlinear fluctuation 
effects - are taken into account). One example of this is 
the cubic symmetry breaking interaction [l^l in the 0{n) 
model, which is relevant by power counting at the Gaus- 
sian fixed point for d < 4, but proves to be irrelevant, 
for sufficiently small n, at the Wilson-Fisher fixed point 
that actually controls the transition for d < 4, at least 
for e = 4 — d sufficiently small. 
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Unfortunately, doing a similar 4 — e analysis of the rele- 
vance of these new nonlinearities in the flocking problem 
would tell us nothing about whether or not these terms 
are relevant in d = 2, since d = 2 is far below the critical 
dimension dc = 4 of the flocking problem. 

Hence, whether or not the exact exponents predicted 
by 01 are correct remains an o pen question. They 
could be; numerical experiments 0,13, [lol, and some real 
experiments [2l| agree with the exponents predicted by 
, which suggests they are, but there is really no way to 
be certain at this point. 

Not all of the predictions of are problematic, how- 
ever. In particular, the claim that long ranged orienta- 
tional order can exist even in d = 2 is unaffected. We 
know this because the nonlinear terms clearly make pos- 
itive contributions to the velocity diffusion "constants" 
^Bcff ^^^'^ ^''^d that they are relevant in the RG sense, 
which means they must change the scaling of the velocity 
fluctuations from that predicted by the linearized theory. 
We know that they are relevant by the following proof 
by contradiction: if all of the nonlinear effects were ir- 
relevant, then simple power counting would suffice to de- 
termine their relevance. But simple power counting says 
that all of the nonlinearities are relevant ioi d < dc ~ 4, 
which contradicts the original assumption that they're 
all irrelevant. Thus, the nonlinearites must change the 
scaling of the velocity fluctuations. Since the effect of 
the nonlinearities is to renormalize the velocity diffusion 
"constants" D^^g and upwards, and since this tends 
to reduce velocity fluctuations, the growth of velocity 
fluctuations with length scale must be suppressed (more 
precisely, its scaZm;; must be suppressed; i.e., it must grow 
like a smaller power of length scale L) than is predicted by 
the linearized version of the equations of motion (|15p and 
p6)) . But those linearized equations predict only log- 
arithmic divergences of velocity fluctuations with length 
scale in d = 2. Hence, the real fluctuations, including 
nonlinear effects, must be smaller than logarithmic by 
some power of length scale, which means they must be 
finite as L — )■ cx). This boundedness of velocity fluc- 
tuations means that long ranged order is possible in a 
two-dimensional flock, in contrast to equilibrium systems 
with continuous symmetries. 

Note that all of the troublesome nonlinearities that 
make it impossible to determine exact exponents in d = 2 
involve the fluctuation Sp of the density p. Therefore, 
if these fluctuations could somehow be "frozen out", it 
would be possible to determine exact exponents in d = 2. 

There are a number of types of flocks in which pre- 
cisely such a freezing out of density fluctuations occurs. 
One class of such systems - flocks with birth and death 
- has been treated elsewhere I idl: others, such as in- 



In conclusion, I have reanalyzed the hydrodynamic the- 
ory of fluid, polar ordered flocks. In addition to iden- 
tifying certain new linear terms in the hydrodynamic 
equations for such systems, which slightly modify the 
anisotropy, but not the scaling, of the damping of sound 
modes in flocks, I have also found that certain nonlin- 
earities that are allowed in equilibrium, and that were 
predicted by earlier work0 to stabilize long ranged or- 
der in d = 2, in fact do not. Other nonlinearities missed 
by earlier work could potentially change the scaling ex- 
ponents from those predicted earlier [7|; whether or not 
they actually do so remains an open question. 
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tality, and to the U.S. National Science Foundation for 
their financial support through awards # EF-lf378f5 
and 1006171. 



compressibl e sy stems |ll|. and systems with long-ranged 
interactions jl2l|. will be addressed in future work[13j. In 
all of these systems, exact scaling exponents can be found 
in d = 2. 
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